Pressure induced superconductor quantum critical point in multi-band systems 
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In multi-band superconductors as inter-metallic systems and heavy fermions, external pressure 
can reduce the critical temperature and eventually destroy superconductivity driving these systems 
to the normal state. In many cases this transition is continuous and is associated with a supercon- 
ducting quantum critical point (SQCP). In this work we study a two-band superconductor in the 
presence of hybridization V. This one-body mixing term is due to the overlap of the different wave- 
functions. It can be tuned by external pressure and turns out as an important control parameter to 
study the phase diagram and the nature of the phase transitions. We use a BCS approximation and 
include both inter and intra-band attractive interactions. For negligible inter-band interactions, as 
hybridization (pressure) increases we find a SQCP separating a superconductor from a normal state 
at a critical value of the hybridization V c . We obtain the behavior of the electronic specific heat 
close to the SQCP and the shape of the critical line as V approaches V c . 



I. INTRODUCTION 

The study of asymmetric superconductivity, i.e., of su- 
perconductivity in systems where different types of quasi- 
particles coexist at a common Fermi surface has raised 
a lot of interest in the last years. This in part is due to 
the relevance of this problem for many different areas in 
physics. It arises in cold atomic systems with supcrfluid 
phases^, in color superconductivity in the core of neutron 
stars^^ and in condensed matter physics^. Furthermore 
it is closely related to inhomogeneous superconductivity, 
as FFLO phases^, since this is a possible ground state 
for asymmetric systems. In condensed matter, as inter- 
metallic materials, due to electrons from different orbitals 
at the Fermi surface, there is a natural mismatch of their 
Fermi wave- vectors. This arises because of the different 
effective masses of these quasi-particles or because they 
occur in distinct numbers per atom. 

Then, in multi-band systems, even in the absence of 
external magnetic fields, one has to consider the possibil- 
ity of inhomogeneous superconductivity or other types of 
exotic ground states as gapless superconducting phases^ 
or phase separation^. 

In this paper we focus on the problem of driving a 
multi-band superconductor to the normal state by ap- 
plying external pressure. When this occurs continuously, 
this transition is associated with a SQCP. The theories 
which have been proposed for the SQCP rely in gen- 
eral on the presence of disorder or magnetic impurities 11 . 
However, there is no reason to expect that this should 
play a role in clean systems driven to a normal state by 
external pressure. In the model we discuss here, mixing 
transfers electrons that participate in Cooper pairing to 
a normal band eventually destroying superconductivity. 
Since there is no dissipation in the normal band, super- 
conductivity disappears due to a loss of coherence in the 
system. Our model considers two hybridized bands in the 
presence of inter and intra-band attractive interactions. 
These interactions are competing and determine the na- 
ture of the zero temperature phase transitions to the nor- 



mal state as hybridization (pressure) increases. We show 
that only when intra-band interactions are dominant this 
transition is continuous. Otherwise it is first order and 
accompanied by phase separation as usual in this case. 

The problem of superconductivity in systems with 
overlapping bands has been treated originally by Suhl, 
Matthias and Walker—. The relevance of the different 
interactions^ has been discussed in terms of an energy 
associated with the Fermi surface mismatch, A(Sk F ) = 
VF^kp and the critical temperature, ksT c . Only when the 
former is much smaller then the latter inter-band interac- 
tions become important. In our approach the mismatch 
5k F depends on hybridization and can be controlled by 
pressure. 



II. MODEL AND FORMALISM 



We consider a model with two types of quasi-particles, 
a and 6, with an attractive inter-band interaction^ g, an 
attractive intra-band interaction U and a hybridization 
term V that mixes different quasi-particles states^. This 
one-body mixing term V is related to the overlap of the 
wave functions and can be tuned by external parameters, 
like pressure, allowing to explore the phase diagram and 
quantum phase transitions of the model. The Hamilto- 
nian is given by 
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where a ka and b ka are creation operators for the light 
a and the heavy b quasi-particles, respectively. The in- 
dex I — a,b. The dispersion relations i L k = fu 

and the ratio between effective masses is taken as a = 



2 



TUa/rrib < 1. For simplicity, we renormalize all the en- 
ergies in the problem by the chemical potential /i a of 
the band a (the non-tilde quantities). Furthermore we 
take h 2 / (2m a Ha) = 1- In this case, the dispersion re- 
lations can be written as, e% = e k /^ a = k 2 — 1 and 

e k = ^fc/Va = "fc 2 - b, With b = [Ibllla- 

The y-term is responsible for the transmutation 
among the quasi-particles. In metallic systems, as 
transition metals 15 , intermetallic compounds and heavy 
fcrmions 16 , it is due to the mixing of the wave- functions 
of the quasi-particles in different orbitals through the 
crystalline potential. In the quark problem, it is the 
weak interaction that allows the transformation between 
up and down- quar k a 2 1 4 ' 1 2 . For a system of cold fermionic 
atoms in an optical lattice with two atomic states (a and 
b), the V-term is due to Raman transitions with an ef- 
fective Rabi frequency which is directly proportional to 
V— . Then, the physical origin of the V-term is differ- 
ent for each case. The main point is that at least in 
inter-metallic systems, hybridization can be easily con- 
trolled by pressure or doping^ allowing to explore their 
phase diagram using these quantities as external param- 
eters. Notice that since hybridization transforms a quasi- 
particle into one another, in its presence only the total 
number of particles is conserved. 

The order parameters that characterize the different 
superconducting phases of the system described by the 
above Hamiltonian are, A ab = —gY^ {b-k-a^ka) and 

A = — UJ2 (b-k-abka) ■ These are related to inter-band 

ka 

and intra-band superconductivity, respectively. The 
anomalous correlation functions can be obtained from 
the Greens functions which also yield the spectrum of 
excitations in the superconducting phases. We use the 
equation of motion method to calculate standard and 
anomalous Greens functions 18 . Excitonic types of cor- 
relations that just renormalize the hybridization^ are 
neglected. The relevant anomalous Greens functions are, 
((a k a;b- k - a )) and {{b ka ;b- k - a )) . When we write the 
equations of motion for them, new Greens functions are 
generated^ 8 -. Some of these are of higher order, as they 
contain a larger number of creation and annihilation op- 
erators than just the two of the initials Greens functions. 
For these, we apply a BCS type of decoupling 18 to reduce 
them to the order of the originals propagators. Finally, 
writing the equations of motion for the new Greens func- 
tions, we obtain a closed system of equations that can 
be solved 5 -. The anomalous propagators from which the 
order parameters are self-consistently obtained are given 
bA 
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where 
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Co = [ejtf (V 2 - A 2 b )]' 



A 2 ef . 



As mentioned before, the poles of these propagators give 
the excitations of the system. Also from the discontinu- 
ity of the Greens functions on the real axis we can ob- 
tain the anomalous correlation functions characterizing 
the superconducting state. In general the appearance of 
exotic superconducting phases is related to the existence 
of soft modes in the spectrum of excitations^. In the 
present case, for the energy of the excitations to vanish, 
it is required that [e£e£ - (V 2 - A 2 ab )} 2 + A 2 ef = 0. 
This can occur by tuning the hybridization parameter, 
such that, V — A ab in which case gapless excitations ap- 
pear at k = k F where = 0. Without this fine tuning 
there are no gapless modes. However, in case the intra- 
band interaction vanishes there is a zero energy mode for 
the wave- vector k, such that, t%e\ — (V 2 — A 2 b ) = 0. We 
will see the effects of this behavior in the next section. 
If, for symmetry reasons, we neglect the term linear in u) 
(Ci =0), we obtain the energy of the excitations in the 
form, 
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The order parameters are determined self-consistently 
by a set of two coupled equations which for finite tem- 
peratures are given by 5 -, 
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where in Bk and uij (fc), we substituted e£ = e and eL 
a + (ae — b). 
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In this equation, we substituted e^, = £ and e£ = 
(e + b — a) I a in and Wj (fc). The quantities p a and 
Pb are the density of states at the Fermi level of the a 
and b bands and lud is an energy cut-off. The right hand 
sides of Eqsl8l and fTOl define the gap functions /(A a 6,A) 
and fb(A a b y A), respectively. In the next section we dis- 
cuss the behavior of the dispersion relations, of the gap 
functions and obtain the free energy of the system. From 
these quantities we obtain the phase diagrams for finite 
and zero temperatures. 



III. RESULTS AND DISCUSSIONS 



Zero temperature phase diagrams for pure inter 
or intra-band interactions 



tion gb(V) is required to sustain superconductivity. At 
this value of the interaction there is a continuous second 
order phase transition and this mixing dependent criti- 
cal interaction #£(14) (or interaction dependent critical 
hybridization) characterizes a superconducting quantum 
critical point. In practice this SQCP can be reached ap- 
plying pressure to the system which is a common proce- 
dure, for example, in the study of HF materials 2 ^. 



B. General case at T ^ 

In this section we consider both intra and inter-band 
interactions and discuss the phase diagrams in the pres- 
ence of controlled mixing and finite temperatures. For 
simplicity, and to show clearly the effect of each term 
we assume strong inter or intra-band terms and in each 
case the other interaction is added perturbatively or ne- 
glected. 

For dominant inter-bands interactions, such that, A — > 
0, we see in figures Q] and [2] that the phase transitions 
that were initially discontinuous, of first order, become 
continuous as temperature increases. 



For completeness we discuss briefly the behavior of the 
system at zero temperature 5 -. For purely inter-band in- 
teractions the transitions are discontinuous and there is 
no SQCP in the system (see Figs. [T]and[2]). From the 
ground state energy we can identify three characteris- 
tic values of the hybridization. Starting from the super- 
conducting ground state, as hybridization increases at a 
value V = V\ appears a minimum in the ground state 
energy at the origin (A = 0) that coexists with an ab- 
solute minimum at finite A a b associated with the super- 
conducting state. Further increasing V there is a first 
order phase transition to the normal state at V = V% , for 
which the energies of the normal and superconducting 
states are degenerate. For still larger V, the supercon- 
ducting state remains as a metastable state until V = V3 
where it stops being a minimum of the energy. The val- 
ues V = Vi, V = V2 and V = V3 for a fixed set of 
parameters (u>d, Pa, 9) yield zero temperature phase di- 
agrams shown in Ref£. It is interesting to point out that 
in the inter-band case sufficiently large values of V can 
give rise to soft modes in the dispersion relations which 
are associated with the presence of Fermi surfaces in the 
superconducting stated. 

For purely intra-bands interactions 2 ^, as in the gen- 
eral case, the dispersion relations of the excitations in 
the superconductor do not vanish for any k or V, since 

the equation [e£e^ — V 2 ] + A 2 e£ 2 = does not present 
any non-trivial solution. For T = , differently from the 
inter-band case, as V increases the value of A the min- 
imum of the ground state energy vanishes continuously 
as the system enters in the normal phase 5 . There are 
no metastable states in this case. We also observe from 
the gap equation that a minimum value of the interac- 
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FIG. 1: (Color online) 3D Graph of the asymmetrical gap 
function /(A a (,, A = 0) versus the inter-band superconduct- 
ing order parameter A a t and temperature. See text for pa- 
rameters used in the graph. 

Figure Q] shows the inter-band gap function as a func- 
tion of the order parameter A a b and temperature. Notice 
that close to where the transition changes from first to 
second order, there is a reentrant behavior which can also 
be seen in figure [3] In figures Q] and we took A = 0, 
V = 0.1, a = 1/7 and b = 0.30. Similar values for the 
parameters were used in the other figures. 

We now treat the case of dominant intra-band inter- 
actions, the inter-band term being considered perturba- 
tively (A a b — > 0) with minor effects. In practice for inter- 
metallic compounds this is the case of greater interest. 
Therefore, we will analyze how the usual intra-band su- 
perconductivity changes under the influence of temper- 
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FIG. 2: (Color online) Graph of the asymmetrical gap func- 
tion versus the order parameter A a b for several values of tem- 
perature and A = 0. These curves are intersections of the 
surface in Figure [1] with planes of constant temperature. The 
parameters are the same used in Fig. [T](see text). 



FIG. 3: (Color online) Renormalized phase diagram showing 
the superconducting critical temperatures as a function of hy- 
bridization, T c x V . The transitions along the critical line 
are continuous and V c is a superconducting quantum critical 
point. The inset shows the behavior of T c close to the SQCP. 



ature and pressure (hybridization). We obtain the vari- 
ation of the electronic term of the specific heat in the 
normal phase as hybridization changes and the system 
goes through the SQCP. 

Figure [3] shows the phase diagram where the critical 
temperature is plotted as a function of hybridization. 
The critical line is a line of second order phase transi- 
tions. For these values of parameters we observe that 
hybridization initially increases the critical temperature 
before destroying superconductivity. Since the transi- 
tions are continuous there is a SQCP at a critical value 
of the hybridization V c . For values of V close to V c , the 
critical line vanishes at the SQCP as T c cx \V — V c y with 
a mean-field shift exponent ip = 1/2 as shown in the inset 
of Fig. El 

In the next section we calculate the linear term of the 
electronic specific heat in the normal phase as the system 
approaches the SQCP. 



of the free energy with respect to the order parameters. 
For illustration we consider the purely inter-band case. 
We get, 



dA ab g + Z^2^ 
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However, sinh (x) / (l + cosh (x)) = tanh (x/2) and mak- 
ing dF/dA ab — 0, we find 
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C. Specific Heat 



and, 



The free energy of the system can be obtained in terms 
of the elementary excitations. It is given by, 
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where anc l w 2 are ^he dispersion relations for A = 
A ab = 0, such that, the function /(w?,^) yields the 
normal contribution to the free energy. 

Notice that instead of obtaining the gap functions from 
the Greens functions they can be found by minimization 



also 



we get, 
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Then, 



where 
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Finally, we get, 
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which is the gap equation obtained previously directly 
from the Greens function. The specific heat is given by, 



Cv (V,T) = -T . 
Using the equation for the free energy, we find, 



(21) 



a k 71=1 



or, 



(23) 

This is shown in figure 2] at the critical value of the hy- 
bridization. 




FIG. 4: (Color online) Electronic specific heat as a function 
of temperature along the critical trajectory (V = V c ). The 
parameters used are given in the text. 

The linearly temperature dependent term of the elec- 
tronic specific heat on the normal phase at very low tem- 
peratures and in particular along the critical trajectory 
is given by, 



7T 2 



(24) 



i=l,2 



where pi (fio) is the density of states of the hybrid bands 
at the Fermi level. This is given by, 
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We wish to calculate pt (ui) at the Fermi surface, i.e., for 
k = kp 1 2 , such that, u>i(k Fl ) = and uj2(kF 2 ) = 0- We 
finally get, 



Pi, 2 (w = Po = 0) 



V 



k Fl2 



2ir 2 



F l,2 



(28) 

The values of k Fl |2 can be easily obtained and we get the 
results for the coefficient j(V) of the linear term of the 
specific heat as a function of hybridization shown in Fig- 
ure [5] We have used the same set of parameters which 
yield the phase diagram shown in Figure [H From these 
figures we notice that the maximum T c occurs in a re- 
gion of the phase diagram where 7 is increasing with V . 
At the critical value of the hybridization the coefficient 
of the linear term is passing through a broad maximum. 
The values of 7 can be measured all along in the nor- 
mal phase, above T c , as hybridization is increasing with 
external pressure and passes through the superconduct- 
ing quantum critical point at V c . A behavior of 7 as 
a function of pressure as that shown in Figure [5] would 
be a strong indication that the present mechanism is re- 
sponsible for destroying superconductivity. Since 7 is 
proportional to the total density of states at the Fermi 
level Figure [5] is helpful to understand the initial increase 
of T c with hybridization as shown in Figure [3] In the 
BCS approximation used here this is related to the den- 
sity of states at the Fermi level and this, as shown by the 
behavior of 7, increases initially as V increases. 

The results of this paper are obtained using a mean 
field approximation which does not include fluctuations. 
The mean-field character of the theory is reflected, for 
example, in the shape of the critical line that vanishes as 
T c oc y/\V - V c \ close to the SQCP as shown in Figure H 
The specific heat calculated above was obtained in the 
normal phase and is due solely to the contributions of 
unpaired quasi-particles in the hybridized bands. There 
are no effects of fluctuations since they are not taken into 
account. 

Due to the nature of the approximations we used, we 
can expect that the superconducting temperatures re- 
flects the variations in the density of states and conse- 
quently in 7 as hybridization is changed. As pointed out 
before the increase in T c for small V is entirely consis- 
tent with the behavior of 7 shown in Figure [5] What 
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FIG. 5: (Color online) Coefficient of the linear temperature 
dependent term of the electronic contribution to the specific 
heat in the normal phase as a function of hybridization and 
at V = Vc for the same parameters as in figure [3] It is also 
shown the value V max for which T c is a maximum in figure 
H For V = V* = V*/fx a = V^t/fJ-a = A one of the hybrid 
bands is above the Fermi level and 7 has a discontinuity. 

is remarkable however is that superconductivity is de- 
stroyed while still both bands contribute to the density 
of states at the Fermi level. It can be easily shown that 

k F 2 ~ fc f\ = 2 \l (^r) 2 + 7T> sucn tnat > hybridization 
increases the mismatch Skp = |&Fi — kp 2 \ between the 
Fermi wave- vectors of the hybrid bands. As It is well 
known^ this can give rise to superconducting instabili- 
ties which however always occur discontinuously through 
first order transitions. In the present case the T = 
superconductor-normal transition is continuous being as- 
sociated with a SQCP. Since there is no dissipation in the 



electronic bands but a lack of coherence in one of them 
we may attribute to this the destruction of the supercon- 
ducting phase. 

It is clear that in inter-metallic systems hybridization 
occurs even at zero pressure. The point we wish to em- 
phasize is that this depends on applied pressure and for 
this reason it can be used as a control parameter. A final 
remark is that as a matter of fact the control parame- 
ter is V/ ix a which is the ratio of hybridization over the 
bandwidth of the light quasi-particles. The latter also 
depends on pressure and although we have naturally as- 
sumed that the ratio increases with pressure, this is not 
immediate and may depend on the particular system. 

We have investigated a mechanism to drive a multi- 
band superconductor to a superconducting quantum crit- 
ical point through the application of pressure or doping 17 . 
It does not rely on the presence of magnetic impurities or 
disorder, but on the sensitivity of hybridization to these 
external parameters. Evidence that the mechanism we 
are proposing is in action can be obtained from mea- 
suring the coefficient of the linear term of the electronic 
specific heat just above the superconducting transition 
as a function of pressure. This shows features which can 
be correlated with the behavior of T c . Superconductivity 
is destroyed when all hybrid bands still contribute to the 
density of states and the zero temperature transition is 
from the superconductor to a metallic state. 
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